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INFALL CAUSTICS IN DARK MATTER HALOS? 
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ABSTRACT 

We show that most particle and subhalo orbits in simulated cosmological cold dark matter halos are 
sm'prisingly regular and periodic: The phase space structure of the outer halo regions shares some of 
the properties of the classical self-similar secondary infall model. Some of the outer branches are clearly 
visible in the radial velocity - radius plane at certain epochs. However, they are severely broadened 
in realistic, triaxial halos with non-radial, clumpy, mass accretion. This prevents the formation of 
high density caustics: Even in the best cases there are only broad, very small (< 10%) enhancements 
in the spherical density profile. Larger fluctuations in p(r) caused by massive satellites are common. 
Infall caustics are therefore too weak to affect lensing or dark matter annihilation experiments. Their 
detection is extremely challenging, as it requires a large number of accurate tracer positions and radial 
velocities in the outer halo. The stellar halo of the Milky Way is probably the only target where this 
could become feasible in the future. 

Subject headings: dark matter — galaxies: formation - Galaxy: halo - methods: n-body simulations 
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1. INTRODUCTION 

Idealized models of self-similar, radial infall onto 
an initial spherical overdensity in an expanding and 
otherwise homogenous univ e rse can be solved exactly 
(|Fillmore fc GoldreichI [l98l iBertschingeij flQSSh (FGB 
hereafter). The two key predictions of these models are: 
(i) Nearly isothermal ( p oc r~'^) halo density p rofiles, as 
found in earlier models ()Gottlll975l : lGunnlll977| i. (ii) Near 
their apocenters mass shells pile up and form spherical 
caustic surfaces of infinite density. This occurs at con- 
stant fractions of the current turnaround radius, result- 
ing in slowly outward moving caustics as the turnaround 
radius increases with time. Each mass shell moves on a 
radial orbit with a nearly fixed (only slowly decreasing) 
apocenter of about 0.85 of its own turnaround radius. 
After a shell falls in it forms part of the first (outermost) 
caustic when the shell is near its first apocenter passage. 
Later, near its second apocenter, it forms part of the 
second caustic, and so on. Such caustics would increase 
the dark matter annihilation r ate (|Bergstrom et al.ll2001l : 
iMohavaee fc ShandarinI 120061 ). and it has been argued 
that infall would lead to peaks in the local CDM velocity 
distribution and change the predict ions for direct dark 
matter detection (jSikivie et al.lll997h . 

Both key predictions of the original FGB-model are at 
odds with the properties of dark matter halos found in 
cosmological simulations: Halo density profiles are found 
to be steeper than p cx r^^ in the outer parts and shal- 
lower in the inner regions (jPubinski fc Carlberd flQQll : 
iNavarro et al.l [T996[ ) . However, by relaxing some of the 
simplifying assumptions of the FGB-model it is possible 
to obtain a more realistic secondary infall model which 
successfully reproduces the halo density profiles found 
in cosmological simulations (e.g. lAscasibar et al.l [20071 
and references therein): allowing for angular momen- 
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tum lowers the inner densities relative to the original 
p oc r^^'^^ profile, simply because fewer particles orbit 
through the central regions (Ryden 1993). More realis- 
tic initial pertu rbations lead to a steeper outer profile 
(jHoffman fc Shaham 1985) . 

Infall does not lead to high density caustic surfaces 
in cosmological simulations. In early simulations most 
particle orbits differed widely from the slowly shrink- 
ing, nearly pe riodic orbits found in the FGB-model 
(jZaroubi et al.lll996) . But these simulations contained 
only a few thousand particles, which leads to signif- 
icantly distorted o rbits due to artificial two body re- 
laxation (Diemand et al.l|2004bfl . Despite the numerical 
noise, IZaroubi et al. ( 1996f ) found that the final particle 
energies are strongly correlated with the initial ones, sug- 
gesting a rather gentle, orderly halo build-up as in FGB, 
and not violent relaxation. Another consequence of this 
fairly regular build-up is that material from biased, early 
high densi ty regions ends up in the inner parts of ha- 
los today (jDiema nd et al.l [2005( 1. However, recent high 
resolution simulations still produce rather smooth mass 
profiles without high de n sity caustics (e. g. iMoore et al.l 
[2(M lHelmi"enil1^00l iDiemand et all l200l . iMoorel 
()2001l ) emphasized that in the hierarchical buildup of 
CDM halos the internal velocity dispersion of infalling 
satellites asat is typically only a few times smaller than 
Chost, which would broaden caustics significantly, to a 
width of about iJhostCsat / Chost- Indeed resolved merg- 
ers do bring in practically all the mass an d there is no 
evidence for smooth accretion (jMadau et a l. 2008). Here 
we show that despite clumpy, non-radial, anisotropic ac- 
cretion into a non-spherical potential, the outer parts of 
halos do show a few weak, broad branches in Vr — r- 
space. Their properties agree well with some of the FGB 
predictions. In contrast to the model however, realistic 
infall " caustics" ^ increase the dark matter densities only 
slightly (~ 10% for the strongest ones). Unfortunately 

* For brevity we refer to the turnaround regions of branches in 
Vr — r-space simply as "caustics", even tough they don't have the 
high densities implied by the true sense of the word. 
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Fig. 1. — Pliase space distribution of subhalos at z=0 separated by the number of orbits tliey liave completed. This figure is also available 
as an mpeg animation running from z=12 to in the electronic edition of ApJ. 



this makes them practically undetectable and irrelevant 
for dark matter detection. 

2. WEAK OUTER CAUSTICS IN Vr - R SPACE 

In this section we compar e the histories o f subh alos in 
the Via Lactea simulation fPieman d et al.l [20071 . VL07 
hereafter)^ with predictions from the FGB-model. Later 
we will include the distribution of all dark matter par- 
ticles. The subhalos are selected by peak circular ve- 
locity Vinax > 3 km s^"'^ before accretion, and are there- 
fore very similar to dark matter particl es in their spa- 
tial a nd velocity distribution (Faltcnbac her fc Diemandl 
120061 ). This sample contains 10,652 halos and subhalos 
with masses above about 3 x 10^ M©. The host halo has 
size of l^nax = 181 km s~^, consistent with the dark halo 
of the Milky Way. The subhalos which have completed 
the same small number of orbits end up in a relatively 
narrow branch in Vr — r space, despite their range of or- 
bital eccentricities ( VL07), the prolate h ost halo poten- 
tial (about 0.8:1, see lKuhlen et al.ll2007D an d occasional 
kicks from larger satellites ([Sales et al ][2007[ ). Each panel 
in Figure [1] would correspond to one infinitesimally thin 
branch in the Vr — r plane of the FGB-model. Already at 
first infall the radial and tangential velocity dispersions 
are quite large: ar — 26 km s~^ and CTtan,2D — 60 km s^^. 
These motions induced by large scale tidal forces and also 
by structures forming within the turnaround region pre- 
vent the formation of thin, high density caustics. 

From the positions of subhalos near their apocenter we 
can measure the position and width of the outer six infall 
caustics (Table [TJ. In the FGB-model the radial velocity 
of a caustic at radius r is Wcaustic(?') = 8/9 r/t, where t 
is the time since the Big Bang^. We define those sub- 

^ Subhalo histories, animated versions of Figures [l] and [2] and 
other movies are available at www.ucolick.org/~diemand/vl 

For simplicity we assume this formula derived by FOB in an 



halos near their k*'' apocenter passage and with a radial 
velocity within 24km s^^ of Wcaustic(?') as the members 
of caustic k. For each member i we then determine its 
turnaround radius tk^i and the ratio of today's position 
to turnaround radius rk^i/tk^i- Note that members of one 
caustic fell in from all directions, not as a bound group, 
and merely share similar turnaround times and radii. 
The scatter in the turnaround radii is relatively large 
and comparable to the width of the caustics. This again 
hints at a broadening of caustics caused already early 
(near the turnaround time) by proper motions. Caustics 
4, 5, and 6 are so wide that they completely overlap, and 
they disappear in the full Vr — r plot of all particles (Fig- 
ure[2]). Caustic 3 is the most prominent one at z=0 and it 
leads to a small over-density of about 10% over a smooth 
density profile around 220 kpc. The abundance of parti- 
cles with \vr — Wcausticl < 24km s~^ increases by a factor 
of about 1.6 around 220 kpc. The outermost caustics 1 
and 2 are quite wide, but parts of their corresponding 
branches from Figure 1 are still visible in Figure [2] Note 
that all of the outer three — r branches have large 
gaps, which are caused by epochs of little or no mass 
accretion. The animated versions of Figures [1] and [2] 
shows how these phase-space features and the gaps in 
between form and move. Whether a certain branch is 
currently visible in the Vr — r plane depends on whether 
there was enough mass falling in at some earlier time to 
reach a sufficient density in the corresponding region of 
phase-space. The clearest caustic at z=0 is the third one; 
it was generated by continuous, large infall of material 
around z=1.7. At later epochs Via Lactea's mass accre- 
tion becomes quite small and sporadic (VL07), leading to 
weak first and second branches with large gaps. The Via 
Lactea halo formed in a series of major mergers before 
z=1.7. These mergers could be responsible for a broad- 

Qrn = 1 universe to be approximately right also in ACDM. 
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TABLE 1 

Positions and widths of caustics and turnaround regions in the Via Lactea halo. 
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0.38 
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0.57 
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0.55 


0.78 


0.52-1.46 


0.832 


0.31 


0.34 
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121 


89-170 


0.67 


157 


105-201 


0.61 


0.81 


0.54-1.46 


0.834 


0.27 


0.30 


Note. 


— Positions and widtlis of tlic outer six caustics and turnaround radii of the 


particles in the caustic tk . Median values 


and 68% ranges 



of all caustic members arc given, A refers to the full width of the 68% range. The number of subhalo members in the caustics 1 to 6 are 551, 300, 
49, 32, 56 and 15. "FGB" refers to predictions from the secondary infall model. 
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Fig. 2. — Phase space distribution of the dark matter particles at z=0. The straight line is I'causticC'')- The lower panel shows the p(r) 
of all particles (black) and of a subset with \vr — ^icausticl < 24km s~^ (red). This figure is also available as an mpeg animation running 
from z=12 to in the electronic edition of ApJ. 



ening of the fourth and higher caustics to a degree where 
they become completely washed out and undetectable in 
the z—0 snapshot. 

3. COMPARISON WITH THE FOB-MODEL 

In Table [T] we make some quantitative comparisons 
with the e = 1 version of the FGB-model, which cor- 
responds to a point-mass as initial over-density: S ~ 
SM/M oc M~^ In CDM the initial peaks are extended. 
The Via Lactea halo has e ~ 0.2 up to about 10^^ Mq. 
For larger masses e increases and it is about one at 
10^^ M0, equivalent to a scale of 150 kpc today. The 
outer halo features discussed here correspond to the e ~ 1 



regime . 

The FGB-model was solved for an Einstein-de Sitter 
{flm = 1) universe, while the Via Lactea simulation as- 
sumed the now standard ACDM model, but we do not 
believe that this difference affects our results: In the 
model, the position of a caustic is similar (0.8 to 0.9) to 
the turnaround radius of the caustic members. Members 
of the first caustic turned around before z=l (earlier for 
higher caustics), when the ACDM model was still matter 
dominated and similar to an Einstein-de Sitter universe. 

^ Larger collapse factors in the inner halo (see Fig. 1 in VL07) 
might be related to the smaller r^/tk ratios in the e = 0.2 FGB- 
model. 
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After turnaround the enclosed mass dominates over the 
enclosed vacuum energy and the dynamics are practically 
the same in both cosmologies. 

We find that the ratios of the current caustic posi- 
tions are in good agreement with the FGB-model (Ta- 
ble [1]). Also the median of the ratios of caustic radii 
to turnaround radii {rk/tk)nicd agree quite well with the 
FGB-model, i.e. the typical orbit s go out clos e to their 
turnaround scale. Earlier models ()Gunn|[T977f ) assumed 
they would decrease by a factor of two during virializa- 
tion. This assumption is also made in the definition 
of the formal virial radius^, which therefore only en- 
closes a frac tion of the material orbiting around a galaxy 
(iPrada et al. 2006, VL07). The {rk/tk)med in the clusters 
from (Diemand et al. 2004aJ) are also close to 0.85, i.e the 
orbits of bound material extend well beyond the formal 
virial radius both in galaxy and cluster halos^. 

The main difference between the FGB-model and cos- 
mological halos is the large scatter in space and veloc- 
ity around the lines in Vr — r space on which all matter 
is found in the model. The scatter prevents the forma- 
tion of high density caustic surfaces and makes the broad 
caustic-like features from cosmological infall practically 
undetectable: To find them one needs an accurate, large 
sample of — r measurements, which are available in 
simulations, but unfortunately not in the real Universe. 
For members of the outer four caustics the scatter in 
the individual rk/tk ratios is consistent with the shape 
of the potential in the outer parts of Via Lactea halo, 
which is about 1:0.8 pro late at z=0 and about 1:0.7 pro- 
late at z=0.5 (jKuhlen e t al. 2007). We also find a few 
outliers with much larger rk/tk ratios, probably resulting 
from dynamical i nteractions between multiple satellites 
(jSales et al.ll2007[ ). The scatter is significantly larger for 
the 5th and 6th caustics, presumably because this mate- 
rial was accreted early (z > 1.7), when a series of major 
mergers caused large fiuctuations in the host potential. 

4. MULTIMODAL RADIAL VELOCITY DISTRIBUTIONS 

From Figure [2] it is clear that the radial velocity dis- 
tribution of particles in a shell in the outer halo will not 
be smooth and unimodal, but rather a superposition of 
modes caused by the various branches in Vr — r space. 
The velocity distribution and relative importance of each 
mode changes with distance from the galactic center, and 
fewer modes are present at large radii (Figure [3]). 

The subhalo radial velocity distributions are similar to 
those of particles and we can use subhalo histories to 
illustrate how the particle distributions divide up into 
their components. Between caustics 2 and 3 one finds 
a broad, doubly peaked distribution. The negative ve- 
locity peak is a superposition of first infall and subha- 
los approaching the center on the second half of their 
first or second orbit. The subhalos on first infall have 
the most negative typical velocities, followed by those on 
their first and second orbit, in qualitative agreement with 
the FGB-model. In contrast to the model, we find these 
components have wide radial velocity distributions which 
largely overlap and give rise to only one broad infalling 

* The virial radius of the Via Lactea halo is rvir = f 104 = 288 
kpc. It encloses a mass of Mvir = 1.5 X 10^2 y^^^ 

® Unlike galaxies halos, typical cluster halos do accrete a lot of 
mass around z = leading to a net in fall between their formal 
virial radius and the turnaround radius llCuesta et al.ll2"007l) . 
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Fig. 3. — Radial velocity distributions of particles in different 
radial shells (solid lines). Only in the inner halo (lower right panel) 
the distribution is uni- modal and close to a Gaussian (dashed line). 
The other panel show multi-modal distributions. To illustrate how 
to decompose these particle distributions we also plot the median 
and 68% ranges of ingoing and outgoing subhalos, which have com- 
pleted 0, 1 or 2 two orbits from bottom to top (circles). On the 
third caustic (lower left panel) the subhalos which have completed 
3 orbits are now near their apocenter. The vertical lines give the 
local caustic velocity f caustic (^) = 8/9 

peak. On the third caustic the broad, double peaked dis- 
tribution from in- and outgoing material is the same as 
above, but now there is an additional peak at low radial 
velocity from particles near their third apocenter. 

In the inner halo particles which completed the same 
number of orbits are distributed so widely that the to- 
tal radial velocity distribution becomes one featureless 
peak. It is close to a Gaussian distribution with zero 
mean velocity and a dispersion of 150 km s^^. But note 
that a smooth, inner f(r,Vr) does not preclude the ex- 
istence of coherent structures in the full phase-space 
density f{x,v), such as bound clum ps (subhalos) and 
their unbound de bris (tidal streams) (jMoore "eral]f200l 
iHelmi et al.l[200l) . 

5. SUMMARY 

We have discovered broad, caustic-like features in 
f{r,Vr) in the outer parts of CDM halos from cosmo- 
logical simulations. They are too broad to cause signifi- 
cant density enhancements and would therefore be very 
challenging to detect. Since a large number of accurate 
distances radial velocities are needed, we believe that the 
stellar halo of the Milky Way is the only system where 
one might possibly detect such features in the foreseeable 
future. 

The basic properties of infall caustics are indepen- 
dent of redshift and halo mass; they are similar in the 
Via Lactea galaxy ha lo and in the cluster halos from 
iDiemand et al.l (|2004af ). The galaxy halos in that work 
were resolved with only a few million particles per halo 
and clearly show similar features. A new, one billion 
particle galaxy halo (Diemand et al. in prep.) also only 
shows broad, weak outer caustics. Since the features de- 
scribed here do not change qualitatively when the nu- 
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mcrical resolution is increased by almost three orders of 
magnitude we conclude that they are well resolved in 
simulations with a few million particles per halo. 
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